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Parametric down-conversion from a wave-equations approach: geometry and absolute 

brightness. 

Morgan W. Mitchell^ 

' ICFO-Institut de Ciencies Fotoniques, Mediterranean Technology Park, 08860 Castelldefels (Barcelona), Spain 

(Dated: 26 September 2008) 

Using the approach of coupled wave equations, we consider spontaneous parametric down- 
conversion (SPDC) in the narrow-band regime and its relationship to classical nonlinear processes 
such as sum-frequency generation. We find simple expressions in terms of mode overlap integrals 
for the absolute pair production rate into single spatial modes, and simple relationships between 
the efficiencies of the classical and quantum processes. The results, obtained with Green function 
techniques, are not specific to any geometry or nonlinear crystal. The theory is applied to both 
degenerate and non-degenerate SPDC. We also find a time-domain expression for the correlation 
function between filtered signal and idler fields. 



I. INTRODUCTION 



Spontaneous parametric down-conversion (SPDC) has 
become a workhorse technique for generation of pho- 
ton pairs and related states in quantum optics. Im- 
provements in both nonhnear materials [l| and down- 
conversion geometries have led to a steady growth in the 
brightness of these sources [3, y, 0, [a, M 0, @, Q ■ 

Apphcations of the bright sources include fundamental 
tests of quantum mechanics, quantum communications, 
quantum information processing, and quantum metrol- 
ogy [10! IHIj Ej [i3| • Although down-conversion sources 
typically have bandwidths of order 10^^ Hz, for the 
brightest sources even the output in a few-MHz window 
can be useful for experiments. This permits a new ap- 
plication, the interaction of down-conversion pairs with 
atoms, ions, or molecules. Indeed, sources for this pur- 
pose have been demonstrated [14| . Many modern ap- 
plications use single-spatial-mode collection, either for 
improved spatial coherence, to take advantage of fiber- 
based technologies, or to separate the source and target 
for experimental convenience. 

Remarkably, despite the importance of bright, single- 
spatial-mode sources, general methods for calculating the 
absolute brightness of such a source are not found in the 
literature. By absolute brightness, we mean the number 
of pairs per second that are collected, for specified beam 
shapes, pump power, filters and crystal characteristics. A 
number of calculations study the dependence of bright- 
ness on parameters such as beam widths or collection 
angles, but these typically give only relative brightness: 
the final results contain an unknown multiplicative con- 
stant [la, [la]- While useful for optimizing a given source, 
they are less helpful when designing new sources. A re- 
cent paper computes the absolute brightness for a specific 
geometry: gaussian beams in the thin-crystal limit [17| . 

In this paper, we calculate the absolute brightness for 



narrow-band, paraxial sources. The results are quite gen- 
eral, for example they apply equally well to crystals with 
spatial or temporal walk-off, for non-gaussian beams, etc. 
The Green-function approach we use is well suited to 
describing the temporal features of the down-conversion 
pairs, and we are able to predict the time correlations in a 
particularly simple way. To our knowledge, this method 
of deriving the time-correlations is also novel. 

Perhaps of greatest practical importance, we derive very 
simple relationships between the efficiency of classical 
parametric processes and their corresponding quantum 
parametric processes. For example, in any given geome- 
try the efficiency of sum-frequency generation and spon- 
taneous parametric down-conversion are proportional. 
This allows the use of existing classical calculations 
and/or experiments with classical nonlinear optics to pre- 
dict the brightness of quantum sources. 

The paper is organized as follows: In section |lTl we de- 
scribe briefly the variety of theoretical treatments that 
have been applied to parametric down-conversion, and 
our reasons for making a new calculation. In Section 
mil we describe the formalism we use, based on an ab- 
stract paraxial wave equation and Green function solu- 
tions. In Section [IV] we calculate the absolute brightness 
and efficiencies for non-degenerate and degenerate para- 
metric down-conversion and the corresponding classical 
processes. In section IVll we summarize the results. 



II. BACKGROUND 



The characteristics of parametric down-conversion light 
have been calculated in a number of different ways. 
Kleinman '18| used a Hamiltonian of the form 



H'^-l /d'xE-xiE 



(1) 



and the Fermi "golden rule" to derive emission rates as a 
function of frequency and angle. Zel'dovich and Klyshko 
[l9l | proposed to use a mode expansion and calculate pair 
rates treating the quantum process as a classical para- 
metric amplifier seeded by vacuum noise . Detailed treat- 
ment along these lines is given in [20|, l2l| . The problem 
of collection into defined spatial modes was not consid- 
ered, indeed the works emphasize that the total rate of 
emission is independent of pump focusing. 

After the observation of SPDC temporal correlations by 
Burnham and Weinberg [22], MoUow [23] described de- 
tectable field correlation functions (coincidence distribu- 
tions) in terms of source-current correlations and Green 
functions of the wave equation. This Heisenberg-picture 
calculation derived absolute brightness for multi-mode 
collection, e.g., for detectors of defined area at defined 
positions. It did not give brightness for single- mode 
collection, nor a connection to classical nonlinear pro- 
cesses. Hong and Mandcl [24| used a mode-expansion to 
compute correlation functions based on the Heisenberg- 
picture evolution and an interaction Hamiltonian of the 
form 



Hr = ^ 



d^x^^^E^E.E^. 



(2) 



As with MoUow's calculation, they find singles and pair 
detection rates, but only for multi-mode detection [30| 
Ghosh, et al. [2^ used the same Hamiltonian in 
a Schrodinger-picture description, truncating the time 
evolution at first order to derive a "two-photon wave- 
function." This last method has become the most popular 
description of SPDC, including work on efficient collec- 
tion into single spatial modes[l^, [la|. Many works along 
these lines are cited in reference [1^. Recently, Ling, et 
al. [l3] calculated the absolute emission rate based on 
a similar interaction Hamiltonian and a gaussian-beam 
mode expansion. In this way, they are able to calcu- 
late the absolute pair rate for non-degenerate SPDC in a 
uniform, thin crystal into gaussian collection modes. As 
described in Sect ion IVCl our calculation agrees with that 
of Ling et al. while also treating other crystal geometries, 
general beam shapes, and degenerate SPDC. 

Notable differences among the calculations include 
Heisenberg vs. Schrodinger picture and calculating in di- 
rect space vs. inverse space via a mode expansion. While 
they are of course equivalent, Heisenberg picture calcu- 
lations are easier to compare to classical optics, while 
Schrodinger picture calculations are more similar to the 
state representations in quantum information. As our 
goal is in part to connect classical and quantum efhcien- 
cies, we use the Heisenberg picture. Also, we note that 
the Schrodinger-picture "two-photon wave-function" has 
a particular pathology: the first-order treatment of time 
evolution means the Schrodinger picture state is not 
normalized and never contains more than two down- 
conversion photons. While this is not a problem for calcu- 
lation of relative brightness or pair distributions [29. |27| . 



it does prevent calculation of absolute brightness. The 
choice of inverse vs. real space calculation is also one 
of convenience: for large angles in birefringent media or 
detection in momentum space, plane waves are the "nat- 
ural" basis for the calculation. However, most bright 
sources use paraxial geometries and collection into de- 
fined spatial modes, e.g., the gaussian modes of optical 
fibers. In these situations, the advantages of a mode 
expansion disappear, while the local nature of the x'^' 
interaction makes real-space more "natural." Thus we 
opt for a real-space calculation. 

Our treatment of SPDC is based on coupled wave equa- 
tions, a standar d ap proach for multi-wave mixing in 
non-linear optics [2^. The calculations are done in the 
Heisenberg picture, so that the evolution of the quan- 
tum fields is exactly parallel to that of the classical fields 
described by nonlinear optics. This allows the re-use of 
well-known classical calculations such as those by Boyd 
and Kleinman [29|. As in the approach of MoUow, we 
use Green functions to describe the propagation, and 
find results that are not specific to any particular crys- 
tal or beam geometry. Unlike MoUow's calculation, we 
work with a paraxial wave equation (PWE). This al- 
lows us to simply relate the classical and quantum pro- 
cesses through momentum-reversal, which takes the form 
of complex conjugation in the PWE. 

We focus on narrow-band parametric down-conversion, 
for which the results are particularly simple. By narrow- 
band, we mean that the bandwidths of the pump and of 
the collected light are much less than the bandwidth of 
the SPDC process, as set by the phase-matching condi- 
tions. This includes recent experiments with very narrow 
filters [ij, but also a common configuration in SPDC, in 
which the down-conversion bandwidth is ^ 10 nm while 
the filter bandwidths are < 1 nm. 



III. FORMALISM 



A. description of propagation 



We are interested in the envelopes £± for forward- 
and backward-directed of parts of the quantum field 
E^^\t,yi) = (£+ exp[-f zA:z] -I- 5_ exp[— ifcz]) exp[— iwi] 
where k is the average wave-number and lo is the car- 
rier frequency. These propagate according to a paraxial 
wave equation 



V±E± = S- 



■±; 



(3) 



where T>± is a differential operator and S± is a source 
term (later due to a x non-linearity). 



The formal (retarded) solution to equation ([3]) is 

£± {x) = £o± [x) + / d^x'g± {x- x')S± {x') (4) 



where x is the four- vector (t,x), £o±{x) is a solution to 
the source-free {S = 0) equation, and Q± are the time- 
forward Green functions, defined by 



V±g±{x;x') ^ 5^{x-x') 
Q±{x;x') =0 t<t'. 



(5) 



For illustration, we consider the paraxial wave equation 
(PWE), for which 



,2 



s± = 



UJ 



9 ' ± 



,(NL) 



(6) 

(7) 



Here V^^ is the transverse Laplacian, k = n(a;)a;/c is the 
wave- number, Vg = duj/dk^ is the group velocity, and 
-p(NL) jg |-]-^g envelope for the nonlinear polarization. 

We note that P± is invariant under translations of x, 
and that time reversal i — > — i is equivalent to direction- 
reversal and complex conjugation, i.e., 'D± —>■ V* . The 
results we obtain will be valid for any equation obeying 
these symmetries. In particular, the results will also ap- 
ply to propagation with dispersion and/or spatial walk- 
off, which can be included by adding other time and/or 
spatial derivatives to 7). 

From the symmetries of 'D±, it follows that the Green 
functions depend only on the difference x — x' , and that 
Q+{t,x;t' ,x') = ^l(i,x', f',x). Also, the time-backward 
(or "advanced") Green functions Ti±, defined by 



V±n±{x;x') = S'^{x-x') 
n±{x;x') =0 t>t' 



(8) 



oheyn±{x;x') ^gi{x',x). 



B. boundary and initial value problems 



If the value of the field is known on a plane z = Zgrc; the 
field downstream of that plane is 

£±{x) = f3z f d^x'g±{x-x')£±{x')5{z' - z,,,) (9) 



where Pz = ±2iA;. Similarly, if the field is known at an 
initial time t — to, the field later is 



£±ix) =Pt I d'^x'g±{x]x')£±{x')6{t' 



-to) (10) 



where Pt = 'iik/vg. Similar relationships hold for the 
advanced Green functions. If the field is known in some 
plane z — zq downstream, then 

£±ix) = PI I d'^x"H±{x]x')£±{x')5{z' - Zq) 

= P* j d^x' £±{x')5{z' ~ za)gUx';x) (11) 

while if the field is known at some time tf in the future, 

£±{x) = P*t I d'^x'n±{x]x')£±{x')5{t' -tf) 

= p; f d^x'£±{x')S{t' -tf)g*Jx';x) (12) 



C. quantization 



The field envelopes are operators which obey the equal- 
time commutation relation 



[£ix,t),£\x',t)]^A^^S^x'~x) 



(13) 



where A-y = y^hu!/2nngeo is a photon units scaling fac- 
tor and Ug = c/vg is the group index. For narrow- 
band fields, A^^ {^^^) describes a photon number den- 
sity, and VgA-^ {£^£) and VgsVgiA~^A-^ (£l£l£i£sj de- 
scribe single and pair fluxes. We find the unequal-time 
commutation relation from equation ([TI 



(14) 



[£{x),£\x')],^,.^PtAl^g{x-x') 



so that {Q\£{x)£\x') |0) = PtA^^g{x]x') for t > t' . For 
the PWE, A-'^Vg = IncEolhuj and PtAi^ = ihjj'^/c^£Q. 

To calculate singles rates, we will need to evaluate ex- 
pressions of the form {££^)- For this, a useful expression 
is derived in the Appendix: Equation (|A2p 



{£ix)£Hx')) 



2hnuj^ 
c^eo 



d'^x"S{z" - Zq) 



xg*{x";x)g{x";x'). (15) 

Here zq is any plane down-stream of x and x' . 

D. single spatial modes 



A single spatial mode M± (x) is a time-independent solu- 
tion to the source-free wave equation X'±M±(x) = 0. 
Af|.(x) is the corresponding momentum- reversed solu- 
tion X';pM|.(x) = 0. We assume the normalization 
J d^x\M±{x)\'^6{z) — 1. For single-mode collection, it 



will be convenient to define the projection of a field £{x) 
onto the mode M as 



(16) 



£M{t) = / d-'xM*{y:)S{z - zo)£ix) 



(here and below, the +/— propagation direction is the 
same for £,M). Here zo is some plane of interest, and 
£m (t) describes the magnitude of the field component in 
this plane. Similarly, if the envelope is constant, the field 
distribution is 

£ix)^£Mit)M{^). (17) 

The optical power is (MKS units) PM{t) — 
2nceo J d^x\£{t,x)\'^d{z - zq) ~ 2nceo\£M{t)\'^- 

Given an upstream source S{x), the M component of the 
generated field is 

£M{t) = I d^xd^x'M*{x)S{z~ zo) 

xg{x;x')Six'). (18) 

If the source is time- independent, then Equation pip and 
the time-translation symmetry of G imply 



£M{t)^j j d?x'M*{^)S{x'). 



(19) 



where g — ~Am{'x)d/c^ , d is the effective nonlinear- 
ity, equal to half the relevant projection of x''^^ ^^id 
Afc = kp — ks — ki is the wave-number mismatch. The 
dimensionless function to(x) describes the distribution 
of X ■ For example in a periodically-poled material 
it alternates between ±1. We can take Afc = with- 
out loss of generality, as the phase oscillation can be in- 
corporated directly in the envelopes. The propagation 
directions (±) will be omitted unless needed for clar- 
ity. Note that for transparent materials x^^-* is real, and 
X^'^\ujp;ujs+uji) = x^^H^s;^p-Wi) = x^^Hwi;^p-^s)- 

First-order perturbation theory is sufficient to describe 
situations in which pairs are produced. For example, if 
£os , £oi , £op are source- free solutions, then 

£s = £os+t^s / d'^x'Qs{x;x') 

xgix')£opix')£lix') + 0{g^). (23) 

and similar expressions for £i , £p are sufficient to give 
the lowest-order contribution to the pair-detection rate 

14^^^) (X (£l£l£i£sj- Higher-order expansions would be 

necessary for double-pair production, etc. 



F. narrow-band frequency Alters 



Similarly, if a product £i (a;i )£^2 {^2 ) is given by a constant 
pair source 5^^^ (a:) as 

fi(a;i)£2(a;2) = / d'^xQi{xi;x)g2{x2;x') 

x5(2)(x'). (20) 

then the time-integrated mode-projected component is 

1 



dti£iMi {tl)£2M2 (^2) 



d^x'M*{^') 

PlzP2z J 

xM;{x')S^^^{x'). (21) 



E. Coupled wave equations 



We now introduce a x nonlinearity, which produces 
a nonlinear polarization that appears as a source term 
in the propagation equations. We consider three fields, 
"signal," "idler" and "pump" with carrier frequencies 
ujs,uJi,ujp and wave- numbers ks,ki,kp, respectively. The 
respective field envelopes £s,£i,£p evolve according to 



UpCp 

Vs£s 
V,£, 



ujpg£s£i exp[iAfcz] 



.2„c- ct.^r^j^fc^i 



ujfg£p£l exp[— iAfcz] 



(22) 



In most down-conversion experiments, some sort of fre- 
quency filter is used. Assuming this filter is linear and 
stationary, the field reaching the detector is 

£'-^Ht) = f dt'F{t - t')£{t') + G{t - t')frcs(t')- (24) 



Here 5res is a reservoir field required to maintain the field 
commutation relations. Assuming the reservoir is in the 
vacuum state, it will not produce detections and can be 

ignored. Defining HpiU.ts) = (£f'\ti)£i'''\ts)) , the 

fields that leave the filter obey 



Hf{U:U) = dt'dt"F,{U-t')Fs{ts-t") 

X {£,{t')£,{t")) . (25) 

In the narrowband case, i.e., when the correlation time 
between signal and idler is much less than the time- 
scale of the impulse response functions, we can take 
{£^{t')£s{t")) « A5{t' - t") where the constant A = 
]dU{£,{U)£s{ts)). We find 

HF{U,ts) « A I dt'F,{U-t')Fs{ts-t') 

= Af{ts-ti). (26) 

With this, we see that the flux of pairs is 

W^'\U-U) = ^^,f'"'"° |.4/(t.-t.)P (27) 



ft2, 



UJsUJi 



with a total coincidence rate of 



A. sum-frequency generation 



W^(2) ^ dUW(^\t,-U) 



2 2 
TlsTliC £"q I . 12 



;i2 






A\' / dt,|2/(i,-t0|^ 



hPujs^i 



\A\^T, 



eff- 



(28) 



We note that 



Tcff 



dr2Ts(r2)Ti(-f7) 



(29) 



where Ts^i(r2) = | J dt exp[iflt]Fs^i{t)\'^ are the signal and 
idler filter transmission spectra, respectively. For this 
reason we refer to Fofr as the effective line-width (in an- 
gular frequency) for the combined filters. Also important 
will be the singles rate 

w^'^ = A;>,,([£f-)(t.)]t£:r^^)(i.)) 

2nsCeo 



dt'dt"F*{ts-t')Fs{ts-t") 
C I dt'\Fs{t,-t')\^ 



2rtsC£o 

risceo 
2fujj, 



CT, 



(30) 



where C = J dt' {£Ut')£s{t")). TeS.s is the effective line- 
width for the signal filter. 



IV. RESULTS 



With the calculational tools described above, we now 
demonstrate the central results of this paper. We first 
express the efficiency of continuous-wave sum-frequency 
generation (SPG) in terms of a mode-overlap integral. 
This effectively reduces the non-linear optical problem 
to three uncoupled propagation problems. We then show 
that the efficiency of parametric down-conversion in the 
same medium is proportional to the SFG efficiency, for 
modes with the same shapes but opposite propagation 
direction. The constant of proportionality is found, al- 
lowing calculations of absolute efficiency based either on 
material properties such as x*"^"* or measured SHG effi- 
ciencies. Similarly, the singles production efficiency is re- 
lated to difference-frequency generation (DFG) and the 
collection efficiency is calculated. The same quantities 
for the degenerate case are also found. 



We consider first the process of SFG, for un-depleted sig- 
nal and idler and no input pump. Signal and idler are 
constant and come from single-modes, 



£Mp{tp) — — i?j\/;£ 



Ms 



C Pz,p 



X / d^x' M;{x')Tn{x')M^{x')Msix') 



— —£]\i,£. 



AluM 



Mi^Ms 



c¥.. 



-ISFG- 



(31) 



The conversion efficiency is 
Pm. 



Q 



SFG 



_ 8L0*npd^\IsFG\' 



2^2 



2w^d 



c'eonpUsrii 



LSFG\ 



(32) 



The efficiency of a cw, single-mode source is thus pro- 
portional to the spatial overlap of the pump, signal, and 
idler modes, weighted by the nonlinear coupling g. 



B. non-degenerate parametric down-conversion 

Next we consider the process of parametric down- 
conversion. Using Equation (P5)l . we can calculate to 
first order in g the correlation function 

{£i{xi)£s{xs)) = ujI j d'^x'Qs{xs,x) 

X (£oAxi)£l,{x'y^ 
xg{x')£o^p{x') 

c^£o J 
xg,{x,,x')g{x')£o,p{x') (33) 

For constant pump and single-mode collection we have 

•^MiMs = j dts {£Mi{ti)£Ms{ts)) 

ihhJs^id 

X I d^x'Ml[x')M*{x')fn{x')Mp{x') 
-£m„Idc- 



ihiUsLUid 



(34) 



We note that Ij^c = ^sfg- Also, the conjugate modes 
describe backward-propagating fields, as if the source 



fields were sent through the nonhnear medium in the op- 
posite direction. Thus if we want to know the brightness 
of down-conversion when all beams are propagating to 
the left, it is sufficient to calculate (or measure) the effi- 
ciency of up-conversion when all beams are propagating 
to the right. Using equations ((32|) and (|34|) we find 



1^ 



Mi Ah 






PpQ 



SFG- 



(35) 



generation. We find 



Pr = 



c^eonsHinp^ ^ 



PpPs / d XiS{zi - zq) 



f3,^, / d^x'g^ix^■,x')mix')Mp{x')M;{x') 



P P 



2luM^ 



r(«) |2 



c-^eonsniUp 
P'pPsQdfg- 



(38) 



C. brightness 



E. singles rates in PDC 



We can now consider the brightness of the filtered, single- 
mode source. The rate of detection of pairs is 



T^(2) ^ 



h?U!sUJi 



^l^preff 



— Fpff— — —PpQ 



H 



P^SFG 



(36) 



This simple expression is the first main result: The rate of 
pairs is simply the joint collection bandwidth Feff, times 
the ratio of frequencies, times the pump power, times the 
up-conversion efficiency QgpQ- Note that the last quan- 
tity can be calculated if the mode shapes and x''^-'(x) 
are known, for example in the paper of Boyd and Klein- 
man, or simulated for more complicated situations. Most 
importantly, it is directly measurable. 



D. difference- frequency generation 



We now consider the classical situation in which pump 
and signal beam are injected into the crystal and idler 
is generated. We will see that this directly measurable 
process is related to the singles generation rate by para- 
metric down-conversion. The generated idler is 

E^{x) = ujf j d'^x'Qi{x]x')g(x')Eop{x')£^Xx'). 

If pump and signal are from modes Mp, Ms, respectively, 
we find 

£i{x) = ^^M^{tp)£*MSts) / d'^x'Qi{x]x') 



xm{x')Mp{x')M*{x'). (37) 



We can find the rate of detection of singles in the mode 
Ms by equation ([50]) and using Equation (|A2p 



C = dts(^£lj^{xs)£Ms{x's] 
= f dUd^x,d^x',Msixs)M:{x',) 

x6{z,s - zo)6{z'^ - zo) {£l{xs)£s{x'^)} 
- I^^^pI''^'- j d^xd^x'Ms[^)g{^)M;{^) 



\Pz 



X (fo.(x)4(^')> M;(x').9(x')Mp(x') 



2hlOiLUgd i2 / i4 lit/ II \ 

n . \£mp\ / dV'(5(z"-zo) 



(3,,, / d3xG,(x";x)M;(x)m(x)Mp(x) 



(39) 



so that 



W^(i) ^ 



CUJiUJs 



32npnsnieo 



Ppl^DFGi ^ cff,j 



— r ff p o^''> 



(40) 



F. conditional efficiency 



The conditional efficiency for the idler (probability of col- 
lecting the idler, given that the signal was collected) is 



1^(2 



Tcff \Isfg\ 



W, 



(1) 



|7-(s) 12 

\-^dfg\ 



(41) 



G. degenerate processes 



The total power generated is Pi = 2cni£o J d^Xi5{zi — 
Zf))\£i{xi)\'^ where zq indicates a plane downstream of the 



Up to this point, we have discussed only non-degenerate 
processes, i.e., those in which the signal and idler fields 



are distinct and do not interfere. This is always the case 
for type-II down-conversion, and will be the case for type- 
I down-conversion if the frequencies and/or directions of 
propagation are significantly different. We now consider 
degenerate processes, in which there is only one down- 
converted field (signal). 

The above discussion is modified only slightly. The signal 
and pump evolve by 

V,£, = io^sg^pSl (42) 



H. second harmonic generation 



The calculation of second-harmonic generation (SHG) 
proceeds exactly as in sum-frequency generation, except 
for the factor of one half and with all "idler" variables 
replaced by "signal" variables. Thus we find 



Pp = P^QSHG 



where 



and 



QsHG 



^P^' \r |2 

A'-shg\ 



2c^eonpn1 



IsHG = /'d3xM;(x)m(x)Af,(x)Af,(x). 



I. average parametric gain 



(43) 



(44) 



(45) 



term is the phase-independent contribution to the gain, 
an experimentally accessible quantity. We have 



= ^^^^^\£0s\'\£p\' J d'xsS{zs - zo) 



Pz,s / d^x'gs{xs;x') 



X m(x')Mp(x')M;(x')| 



c-^n^npEo 
PqsPpQapg- 



-\Iapg\ 



(48) 



J. degenerate PDC 



Next we consider the process of degenerate parametric 
down-conversion, for which 



£s = £as + ^l / d'^x'Qs{x;x') 

^g{x')Eop{x')£l{x'). (49) 



The other classical process of interest is parametric am- 
plification of the signal by the pump. The first-order 
solution for the signal field is 

£s = Eos+ojI I d'^x'Qs{x;x') 

xg{x')Sop{x')£l^{x') 
= ^os+^is. (46) 

The signal power at the output is 

r 

Ps = 2nsceQ / d^XsS{zs ~ ZQ)\£sixs)f 

= 2nsceo / d^XsS{zs - zq) [\£o,s{xs)\'^ 

+2Re[£o^s{xs)£ls{xs)] + \£iAxs)\^) ■ (47) 

The first term is the input signal power Po,s, the second 
term depends on the relative phase (j)p — 2(j)s, and the last 



We find the correlation function 



{£s{xs)£s{x's)) = ojI j d'^x"Qs{xs,x") 

x(£,A<)£Ux"))g{x") 
x£o.p(x") (50) 



at which point it is clear that the only difference from the 
non-degenerate case of Eq. (|55|) will be the replacement 
of idler variables with signal variables. We find 



^^'^ = y.S-r^PpQsHG = -^PpQsHG- (51) 



8 



K. Singles rates (degenerate) 

As before, we can find the rate of detection of singles in 
the mode Ms by equation (PU)) and 

C = I dt',(£{,^{Q£Ms{t's)) 



xS{z'^ - zq)S{z'J - zo) {£l{x',)£s{x'^)) 
X (fo.(x')4(^")) AC(x").9(x")Af,(x") 

X /3,,, f d^xGsix";x)M:ix)m{x)Mp{K) 
The singles rate is thus 

L. Conditional efficiency (degenerate) 



(52) 



(53) 



where m G {s, i,p}, r is the radial component of x, and 
q = z — izR where zr is the Rayleigh range, assumed 
equal for all beams. We assume a periodically-poled ma- 
terial in which x^^^ (z) alternates with period 2t:/Q so and 
we approximate m(x) w eyi^[iQz\dcald. From Equation 
([5T|1 we find 



'SFG 



;- 


TT^ J- 


L/2 
L/2 


-iAkz' 




^y 


2Trr'dr'e^ 


i^ 


, ; 2 




2i 


/ f^pf^sf^i^j^ 




J 


V TT 
r.L/2 




-iAkz' 




^y 


d-' 




^ 




-L/2 iz' 


— i 


Zfl) (i?fcZ' + 


iZR) 



(56) 



where Afc = kp — kg — ki — Q and -Rfc = /j_/fc+ and fc± = 
kp ± (ks + ki). In terms of the dimensionless variables 
K = AkL, C = z'/L, (^ji = Zfi/L we find 



2i 



IsFG = — a/ttAP^A^TT 



where 



T = 



Ci^ 



1/2 



dc 



'IkC, 



27ry_i/2 Mc - *Cfl) (^fcC + <fl) ■ 



(57) 



(58) 



The conditional efhciency is 



W, 



(1) 



-Lcff I^SifGl 



(54) 



From Equation (|32p the upconversion efficiency is then 

n,p Kg Ki 



Q 



SFG 



SttljI 



c^SonpUsrii 



kl 



-Znd^g\T\ 



(59) 



A. Boyd and Klienman, 1968 



EXAMPLE CALCULATIONS 



We now illustrate the preceding, general results with a 
few special cases. We first calculate the overlap integral 
for co-propagating gaussian beams. This allows us to 1) 
compare our results to the classical results of Boyd and 
Kleinman [23|, 2) predict absolute brightness for an im- 
portant geometry, type-II co-linear down-conversion in 
quasi-phase-matched material. Also, we compare to a 
recent calculation of absolute brightness for a specific ge- 
ometry by Ling, et al. |17| . 

We consider collinear, frequency-degenerate type-II PDC 
with circular gaussian beams for signal, idler and pump. 
We take mode shape functions 



M,„(x) 



kmZR 1 

TT q 



^ik^z^ik„,^^ 



(55) 



With this expression we can compare our results to 
those of Boyd and Kleinman [231 f^'' the case of second- 
harmonic generation. As that calculation does not in- 
clude quasi-phase matching, we take Q — 0, and then 
for any reasonable phase-matching we have Up « rig, 
kp « 2ks and thus k- = Rk = 0, k+ w 2kp. We note that 
for i?fc = 0, T becomes equal to the function H of Boyd 
and Kleinman for zero absorption and walk-off angle. We 
find 



Pp = 



Ank^ 



c^eon'sUp 



ZHd'\T\'Pl 



(60) 



Boyd and Kleinman find in Equations (2. 16), (2. 17) and 
(2.20) to (2.24) 



P2^'-^^d^PiLk,^^\Hf 



(61) 



P2 = 






(62) 



When converting this expression to MKS units, cP — 
(P/64TT^eo, and we see that the two calculations agree. 



B. Type-II collinear brightness 



Next we make a numerical calculation for frequency- 
degenerate type-II SPDC, a geometry of current interest 
for generation of entangled pairs, for example. 

The integral T must be evaluated numerically. For a 1 cm 
crystal of PPKTP and a vacuum wavelength A^ = Ai = 
800 nm we have {ns,m,np) = (1.844,1.757,1.964) and 
dcs — 2.4 pm/V so that Rk — 0.04 and the maximum of 
{zr/L)\T\^ « 0.054 occurs at k « -3.0, Ci?, « 0.18. We 
find QsFG = 2.0 X 10~^ W~^ Used as a photon-pair 
source, this same crystal and geometry would yield by 
Equation 



VK(2) = FeffFp- 



:SFG 
16 



(63) 



or a pair generation efficiency of Qsfg/^Q = 0.8 pairs (s 
mW MHz)-^ Note that T^s is the fiher bandwidth in 
angular frequency. 



C. Ling, Lamas-Linares, and Kurtsiefer, 2008 



Recently, Ling et al. calculated the absolute emission 
rate into gaussian modes in the thin-crystal limit of (non- 
periodically-poled) nonlinear material [17| . They arrive 
to a down-conversion spectral brightness of 



dR{uJs) _ / rfcffa^a»£^"^(Afc) \ qj^uj 



(64) 



du)s \ ^ I 27ms rii 

where R is the pair collection rate and 

$(Afc) = I dz f dydxe'^^''Up{r)Us{r)U^{r). (65) 

Here Um describe the mode shapes of the form Um{^) = 
f^ik^z^-{x^+y^)/wl ^^^ Q,^ ^ y^2/nWl are normal- 
ization constants. The field E^ is defined such that 
i?°P — 2apPp/eQnpC where Pp is the pump power, giving 



dRJLJs) 
duj. 



UJs'^i<P 



Trc-^eonpUsni 



Pp \apasa^<^{Ak)\ 



(66) 



Assuming the output is collected with narrow-band fil- 
ters of transmission Tg[u)s),Ti(uJi) for signal and idler, 
respectively, the integrated rate is 



R 



dR{LUs 

diu. 



dnTs{ujp/2 + n)Ti{ujp/2-n) (67) 



where we have assumed dR{ujs)/duJs constant over the 
width of the filters. For comparison, using Eqs ([32|) and 
(1361) . we find 



W-(^) = .^^-^' Pp\Is '^ 



2c^eonpnsni 



FG\ ioff- 



(68) 



Then with Equation (|29|) and noting that in the thin 
crystal limit |/sfg| — \oipCtsOii^{Ak)\, we see that the 
two results are identical. 



VI. CONCLUSIONS 



Using the approach of coupled wave-equations, famil- 
iar from nonlinear optics, we have calculated the abso- 
lute brightness and temporal correlations of spontaneous 
parametric down-conversion in the narrow-band regime. 
The results are obtained with a Green function method 
and are generally valid within the paraxial regime. We 
find that efficiencies of SFG and SPDC can be expressed 
in terms of mode overlap integrals, and are proportional 
for corresponding geometries. Also, we find pair time 
correlations in terms of signal and idler filter impulse re- 
sponse functions. Results for both degenerate and non- 
degenerate SPDC are found. Comparison to classical cal- 
culations by Boyd and Kleinman, and to a recent calcu- 
lation by Ling et al. show the connection to classical 
nonlinear optics and "golden rule" -style brightness cal- 
culations, while considerably generalizing the latter. We 
expect these results to be important both for designing 
SPDC sources, as the results of well-known classical cal- 
culations can be used, and for building and optimizing 
such sources. 



APPENDIX A: ALTERNATE PROPAGATOR 



We can use Equations (|T2l) and p3)l to express the prop- 
agator as 

{£{x)£\x')) = \Pt? I d^x"d^x"'5{t" -tf) 
xSr-tj){£ix")£\x"')) 

xg*{x"-x)g{x"'-x'). 

= \A^Pt\^ j d^x"5{t"-tf) 

y<g*{x"-x)g{x"-x') (Al) 
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Noting that J d-^x"S{t" - tf)g*{x";x)g{x";x') 
Vgjd'^x"6{z" - zo)g*{x";x)g{x";x') we find 



{£{x)£Hx')) 



c^eo 



d''x"5{z" - zo) 



xg*{x"-x)g{x";x') 



(A2) 



APPENDIX B: LORENTZIAN FILTER 



A common filter has a Lorentzian transfer function and 
an exponential impulse response 



terms, not having a filter there at all), filter becomes 



fits - U) 



2 



ti<ts 

exp[-r,(tj-is)/2] ti >ts 



(B5) 



That is, the idler photon will always arrive later, and with 
a distribution (after the signal arrival) that is precisely 
the transfer function of the idler-beam filter. Another 
interesting limit is for matched filters, Fg = F^ = F. 
Then we find 



fits - ti) 



F 



exp[-F|i«-t,|/2]. 



(B6) 



F{t) = -9{t)cM-'Ct/2]. 



(Bl) 



The spectral transmission is Tip.) ^ V'^/iV'^ + 4fi2), i.e., 
unit transmission for constant f , a full- width at half- 
maximum of AilpwHM = r and an area J dflT{Q) = 
TrT/2. If we put a filter of this sort in each arm, the 
output has 



lits-ti) 



r^r. 



dt'eits - t')eiu - 1') 



fir) 



X exp[-F,(t, - t')/2] 
xexphF,(t, -i')/2] 



TsT^ f exp[-F,T/2] t > 
2(F, +F,) lexp[F,T/2] t<0 



(B2) 



The effective bandwidth is 



r.ff^A fdT\f{T)\ 



(B3) 



(B4) 



It is worth noting that in the limit F^ ^ cxd (the limit 
of a broad-band filter in the signal beam, or in practical 



Note that for F^ -^ oo, the detection rate is |^pFi/4, 
i.e., proportional to the idler filter bandwidth F^. The 
reverse, s <-> i is also true, of course. From this we can 
get an idea of the conditional efficiency: The rate for 
filtered signal with any idler is proportional to 



F. > 



F F 



F., 



(B7) 



For example, putting matched filters Fg = F^ = F will 
give a rate proportional to FiFs/(Fs -I- F^) i.e., half of the 
rate without the idler filter. This indicates that, of the 
signal photons that pass the the signal filter, half of their 
"twin" idler photons do not pass the idler filter. 
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